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LONG-TERM  GOALS 


The  proposed  work  investigates  quantitatively  the  interaction  between  wave,  currents  and  seabed  sed¬ 
iments  in  shallow  water  over  a  bed  characterized  by  heterogeneous,  mud-dominated  sediments.  The 
long-term  goal  of  is  to  develop  an  approach  to  characterize  accurately  the  state  of  a  muddy  sea  bed, 
based  on  minimal  prior  information  about  bed  sediment,  and  remote  observations  of  surface  waves  and 
currents. 


OBJECTIVES 

This  work  is  a  collaborative  study  between  University  of  Florida  and  University  of  Delaware  (PI:  T.-J. 
Hsu,  N00014-1 1-1-0272).  The  objective  of  the  project  is  to  investigate  the  possibility  to  predict  bottom 
sediment  processes  using  field  data  collected  during  the  MURI  Wave-Mud  experiment.  The  observa¬ 
tional  data  will  be  used  to  identify  typical,  predictable  scenarios  (sequences  of  states)  of  the  evolution  of 
bed  rheology  under  energetic  waves.  In  parallel,  we  propose  to  develop  a  fully-physical,  high-detail  pilot 
model  for  wave-sediment  interaction.  The  pilot  model  will  have  two  components:  a  numerical  models 
for  small-scale  bottom  sediment  transport  Hsu  et  al.,  2009,  and  a  stochastic  nonlinear  wave  propagation 
(Agnon  and  Sheremet,  1997;  Sheremet  et  al.,  2010).  The  data  will  be  used  to  validate  the  model  and 
propose  simplifications  for  operational  implementations.  These  goals  are  aligned  with  all  three  major 
ONR  research  thrust  areas:  nearshore,  estuarine  and  riverine  processes;  remote  sensing  of  the  coastal 
environment;  and  sediment  transport. 


APPROACH 

Laboratory  and  field  observations  show  that  soft  muddy  bottoms  and  near-bed  fluid  mud  layers  can  dis¬ 
sipate  as  much  as  80%  of  wave  energy  over  a  distance  of  just  a  few  wave  lengths  (Gade,  1957;  Jiang  and 
Mehta,  1995;  deWitt,  1995;  Hill  and  Foda,  1999;  Chan  and  Liu,  2009;  Holland  et  al.,  2009;  and  others). 
Many  theoretical  models  of  wave-mud  interaction  have  been  proposed,  involving  a  range  of  rheologies 
and  dissipation  mechanisms.  Mud  has  been  described  as  a  viscous  Newtonian  fluid  (Dalrymple  and 
Liu,  1978;  Ng,  2000;  deWitt,  1995);  visco-elastic  solid  (Jiang  and  Mehta,  1995);  visco-plastic  Bingham 
material  (Mei  and  Liu,  1987;  Chan  and  Liu,  2009);  or  poro-elastic  material  (Yamamoto  and  Takahashi, 
1985).  Other  processes,  in  addition  to  viscous  dissipation  in  the  mud  layer,  have  been  hypothesized  to 
contribute  to  wave  damping,  such  as  nonlinear  interactions  between  surface  and  interfacial  waves  at  the 
water-mud  separation  surface  (Jamali  et  al.,  2003). 

While  linear  rheology  (simpler,  e.g.,  shear  modulus  or  shear  viscosity  independent  of  strain-rate  ampli¬ 
tude)  is  typically  preferred  in  models  over  complex  of  nonlinear  models  (Chou,  1989;  Chou  et  al.,  1993; 
Mei  and  Liu,  1987),  field  observations  provided  by  the  “Sub-bottom  Field  Experiment”  (N00014-10-1- 
0363)  and  previous  projects  suggest  that  complex  rheologies  might  be  quite  relevant  for  applications. 
A  definite  sequence  of  stages  of  bed  transformation  (see  Section  Figure  5)  emerges  from  the  analysis 
Sheremet  et  al.,  2005;  Jaramillo  et  al.,  2008;  Robillard,  2009;  ?;  ?.  Under  energetic  waves,  the  stiff  bed 
softens,  liquefies,  expands,  and  mixes  with  water.  This  mobilizes  a  surficial  layer  of  sediment  that  is  then 
rapidly  resuspended  by  near-bed  turbulence  (a  burst-like  process),  significantly  increasing  the  suspended 
sediment  concentration  in  the  entire  water  column.  In  turn,  increased  SSC  acts  as  a  negative  feedback 
that  controls  the  further  development  of  the  process  by  dampening  near-bed  turbulence  and  suppressing 
mixing.  As  the  storm  wanes,  decaying  near-bed  turbulence  allows  the  suspended  sediment  to  settle,  lead- 
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ing  to  the  formation  of  fluid-mud  layers.  Eventually,  through  de- watering  and  consolidation,  the  stiff  bed 
state  is  reached  again. 

The  proposed  work  is  based  on  the  following  hypotheses: 

1.  Wave  history  drives  the  evolution  of  the  wave/bed- sediment  system. 

2.  Surface  waves  and  the  muddy  seabed  evolve  as  coupled  systems. 

Wave-induced  bottom  stress  and  turbulence  (that  drive  bed  evolution)  is  balanced  by  mud-induced 
wave  dissipation  (that  controls  bed  stress  levels). 

3.  The  evolution  of  the  bed  follows  predictable  cycles. 

The  wave- sediment  system  is  driven  by  waves,  but  its  evolution  path  is  determined  by  the  coupling 
between  wave  and  sedimentary  processes.  Different  wave  histories  might  therefore  result  in  different 
evolution  paths  for  bed  reworking.  Because  waves  are  the  main  driving  factor,  a  coupled  wave-sediment 
numerical  model  driven  by  wave  evolution  can  be  developed  to  forecast  the  state  of  bed  sediment. 

The  concept  of  typical  scenarios  for  bed-rheology  evolution  during  storms  provides  the  basis  for  inves¬ 
tigating  bed  responses  to  waves  activity.  A  pilot  wave-mud  interaction  model  is  being  build  as  a  loosely 
coupled  wave-mud  interaction  model  based  on  two  existing  codes:  the  numerical  model  for  bottom  sed¬ 
iment  transport  (Hsu  et  al.,  2009),  and  the  stochastic  nonlinear  shallow- water  wave  propagation  (Agnon 
and  Sheremet,  1997;  Kaihatu  et  al.,  2007;  Sheremet  et  al.,  2010).  The  ongoing  work  has  three  directions 
of  research: 

Data  analysis:  reconstruct  the  sequence  of  bed  states  in  storms  captured  in  the  field  observations  and 
identify  typical  scenarios  for  evolution  of  bed  rheology. 

Model  development:  improve  existing  wave  and  sediment  transport  model  and  combine  them  into  a 
coupled  wave-sediment  system. 

Model  validation:  simulate  the  collection  of  storms,  to  investigate  the  limitations  of  the  model,  as  well 
as  to  understand  typical  evolution  patterns  of  bed  rheology. 


WORK  COMPLETED 

Field  experiment  and  data  analysis:  The  “Sub-bottom  Field  Experiment”  project  provided  informa¬ 
tion  about  the  evolution  of  the  bed  under  wave  action.  In  support  the  data  collection  effort  of  the 
Traykov ski/Trowbridge  group,  the  project  deployed  instrumentation  to  capable  of  high-resolution  mea¬ 
surements  of  full  water  column  hydrodynamics  and  near-bed  sediment  dynamics.  An  example  of  an 
instrumented  platform  is  shown  in  Figure  (Figure  2).  To  investigate  the  sediment  mass  exchange  be¬ 
tween  the  bed  and  the  water  column,  a  new  method  for  estimating  the  vertical  profile  of  suspended 
sediment  concentration  was  developed  and  tested  using  field  observations  collected  in  2008  and  2010 
(Sahin  et  al.,  2012b). 

Model  development:  Both  numerical  components  of  the  pilot  numerical  model  -  stochastic  nonlinear 
wave  propagation  (Agnon  and  Sheremet,  1997),  and  small-scale  sediment  transport  (Hsu  et  al.,  2009), 
were  implemented  and  thoroughly  tested  on  field  observations  collected  during  the  2008  and  2010  field 
experiments  (e.g.,  Safak  et  al.,  2010;  Sheremet  et  al.,  2010;  Sahin  et  al.,  2012a,b).  A  typical  bed  evolution 
cycle  has  been  identified  and  studied  in  detail  (e.g.,  Sahin  et  al.,  2012b,  see  next  section).  An  analysis  of 
bed  reworking  during  all  observed  storms  is  ongoing. 


2 


When  the  fluid  mud  model  was  originally  developed  Hsu  et  al.  (2009),  there  was  no  available  information 
on  the  mud  rheological  properties  and  an  empirical  closure  used  in  Le  Hir  et  al.  (2001)  was  adopted. 
Consequently,  the  rheology  closure  used  in  Hsu  et  al.  (2009)  only  allows  qualitative  understanding  on 
the  effect  of  rheology  on  fluid  mud  transport.  Without  detailed  the  information  on  mud  rheology,  it  is 
not  possible  to  rigorously  validate  the  numerical  model  and  to  develop  predictive  capability.  Hsu  and  his 
visiting  scholar  Mr.  Wen- Yang  Hsu,  a  PhD  student  of  National  Cheng-Kung  University,  Taiwan,  further 
investigated  this  issue  (see  below).  Mr.  Wen- Yang  Hsu  was  awarded  a  prestigious  9-month  scholarship 
from  National  Science  Council  of  Taiwan.The  new  formulation  on  rheological  closure  will  be  used  in 
the  numerical  model  to  simulate  field  experiment. 

RESULTS 

Field  Experiments:  Figure  3  shows  and  example  of  observations  of  waves  and  bed  state  response  (this 
case:  a  frontal  storm  observed  in  March  2008).  The  event  is  associated  with  fairly  energetic  southward 
winds  and  currents  which  seem  to  be  due  to  a  superposition  of  low  tide  and  the  flushing  of  the  coastal 
setup  post-frontal  storm  passage.  The  sea-floor  response  can  be  inferred  from  the  PC-ADP  acoustic 
backscatter  (Figures  3c),  based  on  the  location  of  maximum  intensity.  At  P3,  bed  elevation  changes 
by  10-15  cm,  consistent  with  previous  observations  (Jaramillo  et  al.,  2008;  Sheremet  et  al.,  2010),  and 
suggesting  a  significant  bed  reworking  and  wave-sediment  coupling. 


Wave  dissipation  mechanisms:  An  inverse  modeling  approach  based  on  the  nonlinear  wave  model 
developed  by  Agnon  and  Sheremet  (1997)  was  used  to  investigate  net  wave  dissipation  (Sheremet  et 
al.,  2010,  Figure  4).  The  dominant  wave  dissipation  mechanism  is  wave-bottom  interaction,  and  that 
the  process  is  triggered  by  the  reworking  of  bed  sediment  by  waves.  Wave  dissipation  increases  during 
a  storm,  as  the  bed  is  softened  by  waves  and  sediment  is  re-suspended.  Maximum  of  mud-induced 
dissipation  (about  50%  loss  of  incoming  energy  flux  loss  over  approx.  4-km  propagation  distance)  is 
observed  at  the  end  of  the  storm,  over  under-consolidated  (gelling)  bed  sediment  (Figure  5b,  8a). 

The  contribution  of  the  nonlinear  interactions  is  expressed  in  transfers  from  the  peak  of  the  spectrum 
toward  higher  and  lower  frequencies,  resulting  in  a  increased  apparent  dissipation  of  the  spectral  peak 
and  net  growth  in  the  high  and  low  frequency  bands.  This  trend  is  not  captured  by  the  linear  model,  which 
suggests  that  neglecting  the  effect  of  wave  nonlinearity  can  lead  to  aliasing  nonlinear  energy  transfers 
into  dissipation  effects,  distorting  the  representation  of  mud-induced  dissipation. 


Bed  Evolution:  Under  energetic  waves  (Figure  5,  P3)  the  stiff  bed  softens,  liquefies,  expands,  and 
mixes  with  water.  The  mobilized  surficial  layer  of  sediment  is  then  rapidly  resuspended  significantly 
increasing  the  suspended  sediment  concentration  in  the  entire  water  column.  Increased  (SSC-induced) 
stratification  controls  the  process  by  dampening  near-bed  turbulence  and  suppressing  mixing  (Safak  et 
al.,  2010;  Sahin  et  al.,  2012a).  As  the  storm  wanes,  decreasing  near-bed  turbulence  allows  settling, 
and  the  formation  of  fluid-mud  layers  (Safak  et  al.,  2010;  Sheremet  et  al.,  2010).  Eventually,  through 
de-watering  and  consolidation,  the  stiff  bed  state  is  reached  again. 


Vertical  profile  of  suspended  sediment  concentration:  The  new  method  was  proposed  by  Sahin  et  al. 
(2012a)  for  estimating  the  vertical  profile  of  SSC  based  on  the  acoustic  backscatter.  The  method  was  used 
to  reconstruct  the  evolution  of  the  suspended  sediment  mass  (Figure  6b).  In  addition  to  flow  velocity, 
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the  SSC  profiles  can  be  used  to  constrain  the  fine-scale  sediment  transport  model  of  Hsu  et  al.  (2009). 
In  turn,  the  model  output  provides  quantities  that  are  difficult  to  observe  directly,  such  as  the  near-bed 
turbulent  stress  (Figure  6d).  The  numerical  model  suggests  a  threshold  bed-stress  value  of  0.52-0.75  Pa 
for  bed  mobilization  (Atchafalaya  mud).  This  value  is  itself  dependent  of  the  wave  activity  history. 

Bed  hysteresis  and  mass  balance:  The  storm  cycle  is  analyzed  in  Figures  7  and  8  (?)  in  terms  of 
the  mass  exchange  between  bed  and  the  observable  water  column  (the  first  approximately  1  mab).  The 
transfer  function  that  maps  the  the  water-column  mass  flux  (rate  of  suspended  sediment  mass  change  per 
unit  time  and  unit  horizontal  area)  onto  the  bed  volume  flux  (rate  of  volume  change  per  unit  time  and  unit 
horizontal  area)  is  the  bed  concentration.  If  during  the  event  the  bed  concentration  were  constant  and  the 
evolution  of  the  bed  and  water  column  sediment  content  were  controlled  strictly  by  bed/water-column 
mass  exchange,  the  shapes  of  the  two  curves  would  be  identical. 

An  estimate  of  the  bed  density  derived  by  directly  dividing  the  two  fluxes  is  shown  in  Figure  8  (circles). 
The  density  has  singularities  for  the  zeros  of  the  bed  flux  (i.e.,  corresponding  to  moments  when  the  water- 
column  mass  changes  but  the  bed  position  is  constant).  Because  the  estimate  is  also  distorted  by  possible 
lateral  advection  of  sediment,  the  periods  with  significant  advective  contribution  (green  segments,  Figure 
7  and  8)  should  be  discarded.  (Advection  however,  is  estimated  to  contribute  less  than  10%  of  the 
suspended  sediment  mass.)  During  the  event,  the  density  of  the  bed  stays  consistently  lower  that  the 
gelling  density  for  Atchafalaya  mud,  but  within  ranges  that  cannot  be  described  as  a  purely  Newtonian 
viscous  fluid.  Toward  the  end  of  the  event,  the  density  increases  steadily  toward  and  past  the  gelling 
point,  consistent  with  inverse  modeling  results  reported  by  Sheremet  et  al.  (2010). 

Modeling  of  Wave-bed  coupling:  Field  observations  of  the  vertical  structure  of  the  flow  velocity  during 
the  March  2008  storm  suggest  that  a  surficial  bed  layer  as  thick  as  20  to  30  cm  oscillates  with  the  waves 
(maximum  RMS  horizontal  velocity  observed  during  the  most  active  period  is  of  the  order  of  20  cm/s) 
while  also  sliding  down-slope  with  a  mean  velocity  of  the  order  of  5  cm/s  (similar  to  that  observed  by 
Jaramillo  et  al.,  2008).  The  fact  that  wave-induced  oscillations  of  the  bed  (which  in  turn  contribute  to 
wave  dissipation)  are  associated  with  high-density  mud  layers  suggests  that  complex,  nonlinear  mud 
rheologies  (visco-elastic/plastic  or  Bingham  type)  might  play  a  role  in  modeling  wave-bed  interaction 
(mud-induced  wave  dissipation,  and  wave  -induced  bed-reworking). 

Most  numerical  model  of  wave-mud  interaction  lack  of  detailed  model-data  comparison,  especially  con¬ 
current  validation  of  free-surface  damping  rate  and  velocity  profiles  in  the  mud  layer.  As  part  of  his  PhD 
study,  Mr.  Wen- Yang  Hsu  has  carried  out  a  comprehensive  series  of  laboratory  experiments  on  wave 
propagation  over  a  bottom  mud  layer  (kaoline  mud).  In  addition  to  the  surface  wave  field  and  bottom  ve¬ 
locity  profiles  in  the  mud  layer,  detailed  rheological  properties  (viscosity  and  yield  stress)  of  the  kaoline 
mud  are  also  measured  via  a  rheometer.  Measured  rheological  properties  were  used  in  the  2DV  nu¬ 
merical  model  of  wave-mud  interaction  (Torres-Freyermuth  and  Hsu  2010),  which  solve  mud  transport 
and  nonlinear  wave  propagation,  to  simulate  the  flume  experiments  and  demonstrate  good  agreement 
between  the  model-predicted  and  measured  wave  shapes  and  dissipation.  Experimental  data,  rheology 
data,  linear  model  of  Dalrymple  and  Liu  (1978),  and  the  nonlinear  model  of  Torres-Freyermuth  and  Hsu 
(2010)  were  adopted  to  carry  out  forward  and  backward  modeling  of  wave-mud  interaction.  Measured 
rheology  of  kaoline  exhibits  hybrid  properties  of  Bingham  and  pseudo-plastic  fluid.  Moreover,  mea¬ 
sured  time-dependent  velocity  profiles  in  the  mud  layer  revealed  that  the  shear  rate  under  wave  loading 
is  highly  phase-dependent.  Measured  shear  rate  and  rheological  data  allow  for  back-calculating  the  time- 
dependent  viscosity  of  the  mud  layer  under  various  wave  loading,  which  also  fluctuates  up  to  one  order 
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of  magnitude  during  one  wave  period.  However,  the  resulting  time-dependent  bottom  stress  was  shown 
to  fluctuate  only  within  25%  of  its  mean.  Measured  wave-averaged  bottom  stress  was  correlated  with 
wave  damping  rate  in  the  intermediate  wave  energy  condition. 

Commonly  adopted  constant  viscosity  assumption  was  then  evaluated  via  linear  and  nonlinear  wave- 
mud  interaction  models.  When  driving  the  models  with  measured  wave-averaged  mud  viscosity  (forward 
modeling),  wave  damping  rate  was  generally  over-predicted  for  low  wave  energy  condition.  On  the  other 
hand,  when  a  constant  viscosity  is  chosen  to  match  the  observed  wave  damping  rate  (backward  model¬ 
ing),  the  predicted  velocity  profiles  in  the  mud  layer  are  not  satisfactory  and  the  corresponding  viscosity 
is  lower  than  the  measured  value.  These  discrepancies  are  less  pronounced  when  waves  become  more 
energetic.  Differences  between  the  linear  and  nonlinear  model  results  become  significant  in  low  energy 
condition,  suggesting  an  amplification  of  wave  nonlinear  effect  due  to  non-Newtonian  rheology.  In  gen¬ 
eral,  the  constant  viscosity  assumption  for  modeling  wave-mud  interaction  is  only  appropriate  for  more 
energetic  wave  condition.  A  manuscript  summarizes  our  main  finding  was  submitted  for  publication. 


Bed  Sediment  Processes:  ID  vs.  3D  Modeling  Observations  and  numerical  simulations  (Sheremet  et 
al.,  2010)  show  that  the  various  transient  states  of  the  bed  during  this  reworking  cycle  are  relevant  for 
wave  dissipation  processes  (Dalrymple  and  Liu,  1978;  Mei  and  Liu,  1987;  Maa  and  Mehta,  1990;  Jiang 
and  Mehta,  1995) 

The  process  of  sediment  mobilization  through  bed  reworking  by  hydrodynamics  suggests  that  the  seabed 
might  play  the  role  as  a  sediment  source.  This  alternative  was  investigated  by  Safak  et  al.  (2010);  Sahin 
et  al.  (2012a,b),  using  field  observations  collected  during  the  ONR-supported  MURI  wave-mud  experi¬ 
ments  of  2008  and  2010.  In  their  preliminary  analysis,  Sahin  et  al.  (2012b)  used  the  ID  vertical  model 
developed  by  Hsu  et  al.  (2009).  They  note  that  the  sediment  balance  is  typically  dominated  by  the  mass 
exchange  between  the  bed  and  the  water  column,  with  the  important  exception  of  water  fronts  passing 
over  the  observation  site,  likely  associated  with  of  freshwater  plumes  originating  from  the  nearby  river. 
The  fronts  produce  a  transient  lateral  (horizontal)  convergence  of  sediment-flux  (Figure  ),  that  cannot  be 
captured  by  a  ID  vertical  model.  Sahin  et  al.  (2012b)  hypothesize  that  water  fronts  are  shore-duration 
processes,  and  overall  the  balance  remains  ID,  even  when  the  plume  passes  over  the  site. 

The  apparent  contradiction  between  the  good  performance  of  ID  models  and  the  existence  of  strong-but- 
transient  3D  processes  such  as  fronts  pose  an  important  problem  for  sediment  transport:  What  are  the 
dominant  sources  of  sediment  on  a  shallow  shelf,  and  under  what  conditions  are  they  active?  The  issue  is 
important  for  modeling  and  prediction,  obviously,  because  ID  modeling  is  lighter/faster,  and  local  -  i.e., 
requires  only  knowledge  of  the  local  flow/sediment/bed  conditions.  However,  ID  modeling  is  predicated 
on  the  assumption  of  non-essential  horizontal  flux  convergence. 


IMPACT/APPLICATIONS 

Much  of  the  present  and  near-future  Navy  capability  on  predicting  regional  and  nearshore  processes 
assumes  a  sandy  (non-cohesive)  sedimentary  environment.  The  present  research  enhances  this  capabil¬ 
ity  by  providing  field  data  essential  for  model  validations  and  by  identifying  processes  and  developing 
mechanisms  which  allow  expansion  into  areas  with  significantly  different  characteristics. 

One  of  the  direct  implications  of  the  present  research  is  the  developing  the  foundation  for  the  develop¬ 
ment  of  a  coupled  hydrodynamic-seafloor  prediction  model  for  muddy  environments. 
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RELATED  PROJECTS 


The  project  represents  a  convergence  of  several  directions  of  research  (near-shore  wave  modeling,  cohe¬ 
sive  sediment  transport,  the  development  of  operational  forecasting  tools  for  near-shore  circulation  and 
waves,  increase  use  of  remote  sensed  information,  etc)  and  etc),  and  collaboration  efforts  circumscribed 
by  the  MURI-lead  effort  to  understand  wave-mud  interaction. 

The  field  experiment  is  coordinated  in  collaboration  with  other  MURI  related  projects.  The  scope  and 
approach  of  the  present  research  builds  on  the  strong,  ongoing  collaboration  between  U.  Florida  and  U. 
Texas  and  U.  Delaware,  illustrated  by  a  number  of  papers  in  print  and  in  preparation.  The  field  work  was 
coordinated  with  with  the  MURI  group  of  researchers,  especially  regarding  observational  data  sharing 
(boundary  layer  and  sediment  characteristics,  Traykovski,  Kineke,  Dalrymple),  and  other  researchers  that 
participated  in  the  MURI-lead  field  experiment  (Elgar,  Raubenheimer,  Allison).  The  work  represents  a 
natural  continuation  and  expansion  of  the  Pis  ongoing  research  projects.  The  proposed  work  also  builds 
on  our  previous  collaborations  on  wave  modeling  with  Kaihatu  (Texas  AM). 

This  research  also  benefits  from,  and  enhances,  parallel  research  (Sheremet)  funded  under  NOPP  to 
improve  existing  operational  wave-forecasting  systems  (Wave Watch  III,  SWAN,  etc)  by  developing  and 
implementing  numerical  modules  for  wave-mud  interaction  and  nonlinear  waves  physics. 

The  bottom  boundary  layer  fluid  mud  modeling  component  of  the  proposed  work  also  benefits  from, 
and  enhances  parallel  research  (Hsu)  funded  by  ONR  to  develop  multidimensional,  turbulence  resolving 
model  for  fine  sediment  transport  driven  by  waves  and  tidal  currents. 
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Figure  1:  a)  Plan  view  of  the  Atchafalaya  shelf  showing  the  location  of  the  experiments  conducted  in 
2008  (light  red  crosses,  platforms  Pl-3)  and  2010  (MURI)  red  cross,  b)  Magnified  area  of  the  2010 
MURI  experiment  with  the  locations  of  the  three  MURI  platforms  ( blue  triangles )  and  Sheremet  & 
Allison  array  ( red  crosses).  An  ADCP  and  a  pressure  sensor  were  deployed  farther  offshore  (approx. 
18-m  depth )  to  provide  boundary  conditions  for  wave  propagation. 
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Figure  2\  Left:  An  instrumented  platform.  Instrumentation  includes  a  PC-ADP  (A),  an  ADCP  (B), 
an  ABS  (C),  a  CT probe  (D),  turbidity  sensors  one  OBS-5  (F),  and  two  OBS-3  (E,  partially  visible 
behind  the  OBS-5),  and  an  acoustic  pinger  (G).  Right:  Pore-pressure  array  ready  to  be  deployed 
(Spring  2010)  and  his  designer,  Uriah  Gravois  (U.  Florida  graduate  student).  Black  cylinder 
contains  electronics  ( Onset  Computer  Corporation  Tattletale  8  Data  logger,  Persistor  Memory 
Expansion  (Paroscientific  Pressure  Sensor  included  in  the  housing).  The  white  cylinder  is  the  probe, 
(4 pore-pressure/thermistors).  Two  Sontek  Hydra  ADVs  and  their  battery  canister  (large  white 

cylinder)  can  also  be  seen. 
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(a)  PI  (29  deg  11.815,  91  deg  36.731  W),  7-m  depth.  (b)  P3  (29  deg  15.574,  91  deg  34.267W),  4-m  depth. 


Figure  3  :  Observations  of  waves  and  suspended  sediment  concentration  (SSC)  at  two  platforms  (PI 
and  P3,  see  Figure  1).  (a)  Significant  wave  height  of  sea  (blue,  f  >0.2  Hz)  and  swell  (red,  f  <0.2  Hz) 
bands.  Multi-color  curve  shows  the  wind  speed  and  direction.  ( b )  Normalized  spectral  density  of  the 
sea  surface  elevation,  (c)  Peak  wave  propagation  direction  for  each  frequency  band  in  the  power 
spectrum  (for  both  winds  and  waves,  the  directions  indicate  where  the  flow  is  toward,  i.e.,  N  means 
toward  North).  The  wave  directions  are  shown  only  for  frequencies  with  spectral  density  above  some 
“ significance  threshold”  (arbitrary),  (d)  Normalized  acoustic  backscatter  records  of  the 

downward-looking  PC -ADP. 
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Spectral  density 


Figure  4.  Example  of  numerical  simulations  of  wave  spectrum  evolution  (a,c)  and  dissipation  rates 
(b,d)for  the  storm  of  March  10-llth  2006  (red  -  observations;  blue  -  model;  black  dashed  - 
mud-induced  dissipation  rate,  Ng,  2000;  crosses  -  net  “ linear ”  dissipation  rate,  including  wind  input, 
white-capping,  and  mud-induced  dissipation).  The  nonlinear  transfer  of  energy  from  the  peak 
toward  higher  and  lower  frequencies  appears  to  increase  the  net  dissipation  of  the  spectral  peak  and 
results  in  net  growth  rates  for  higher  and  lower  frequency  bands.  Nonlinear  wave-wave  interaction 
conserves  energy,  therefore  it  does  not  contribute  to  the  bulk  (frequency  integrated)  wave 

dissipation/growth. 
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Figure  5:  Analysis  of  PC-ADP  backscatter  showing  a  20-30-cm  thick  surficial  layer  of  the  bed 
oscillating  and  sliding  down-slope  during  the  storm,  a)  Significant  wave  height  (blue:  short  waves, 
red:  swell);  and  wind  speed  and  direction  ( color  code  indicated  the  direction  the  wind  blows  toward), 
b)  Normalized  PC-ADP  backscatter  intensity.  The  lines  represent  the  location  of:  maximum 
backscatter  intensity  ( triangles );  zero  mean  horizontal  velocity  (stars),  and  zero  RMS  horizontal 
velocity  (circles).  The  continuous  thick  line  is  a  smoothed  estimate  of  the  bed  position.  A  surficial 
bed  layer  of  approximately  20-30  cm  oscillates  with  the  waves  and  slides  down-slope.  Arrows  mark 
the  hypothesized  stages  of  bed  evolution:  (A)  solid  bed;  ( B )  breaking  of  the  bed  matrix  and  water 
absorption  (liquefaction/fluidization/expansion);  (C)  bed  erosion;  (D)  settling  and  bed  accretion;  ( E ) 
formation  of  fluid  muds.  The  process  is  followed  by  eventual  de-watering/consolidation  (not  shown). 
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Figure  6:  A  reconstruction  of  the  evolution  of  the  vertical  profile  of  the  suspended  sediment 
concentration  during  the  storm  of  March  2008.  a)  estimates  based  on  inverting  the  PC-ADP 
backscatter  intensity  (Sahin  etal.,  2012a),  and  numerical  simulations  using  the  model  of  Hsu  et  al. 
(2009)  ofb)  suspended  sediment  concentration  profile,  and  c)  bottom  shear  stress. 
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Figure  7:  Sediment  mass  balance  between  the  water-column  and  the  sea  bed  during  the 
bed-reworking  cycle:  sediment  flux  vs.  bed  stress,  a)  water  column:  ter-column  and  bed  sediment 
flux  vs.  bed  stress.  The  colors  marks  time  intervals  that  can  be  identified  as  dominated  by 
bed/water-column  mass  exchange,  or  by  lateral  advection.  The  transfer  function  that  maps  the  curves 
in  panels  a)  and  b)  into  each  other  is  the  bed  concentration.  If  during  the  event  the  bed  concentration 
were  constant  and  the  evolution  of  the  bed  and  water  column  sediment  content  were  controlled 
strictly  by  bed/water-column  mass  exchange,  the  shapes  of  the  two  curves  would  be  identical. 


1300 


1250 


E  1200 

O) 


%  1150 
c 
5 
T3 


•o 

9 

CO 


1050 


1000 


|  Adv.Gain 
|  Bed  Gain 
I  Bed  Loss 


(a) 


E 

cn 


0.8 
0.6 
0.4 
0.2 

0  _ 
X 

-0.2  g 
o 

-0  4  w 
V) 

-0.6 

-0.8 

-1 


03/18  04/00  04/06  04/12  04/18  05/00  05/06 
Time  (day/hour) 

Figure  8:  a)  Estimates  of  bed  concentration  vs.  time,  b)  Corresponding  water-column  hysteresis 
curves  (compare  with  Figure  7a).  The  circles  represent  the  bed  concentration  values  derived  directly 
from  Figure  7,  by  dividing  the  two  fluxes;  lines  represent  estimates  based  on  various  degrees  of 
smoothing  of  the  direct  estimate.  In  a),  the  dashed  line  marks  the  gelling  density  for  the  Atchafalaya 

mud  (Robillard,  2009). 
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Figure  9:  Laboratory  experiment  of  wave-mud  interaction  (in  collaboration  with  W.-Y.  Hsu  and 
H.-H.  Hwung,  NCKU,  Taiwan )  reveal  insights  in  the  mud  layer  driven  by  surface  waves  (depth=24 
cm,  wave  height=4  cm,  mud  thickness  =6  cm).  (A)  Measured  velocity  profile  during  flow  reversal 
agrees  with  nonlinear  model  results  ( red  curve).  However,  model  results  are  sensitive  to  the  viscosity 
used  ( see  blue  curve  when  viscosity  is  increased  by  5  times).  (B)  Measured  viscosity  in  the  mud  layer 
(see  B2)  changes  by  one  order  of  magnitude  during  the  wave  passage  (see  (Bl)  for  free-stream 
velocity).  However,  measured  shear  stress  only  changes  within  20%  of  its  mean  value  (see  B3). 
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Figure  10:  Observations  of  the  passage  of  a  fresh-water  front  over  the  experiment  site  during  the 
storm  of  March  3rd  to  5th,  2008.  a)  Wind  speed  and  direction,  and  significant  wave  heights 
( short-wave :  blue  and  swell:  red);  b)  Salinity  (blue)  and  temperature  (red)  at  55  cmab;  c)  Vertical 
profile  of  S SC  estimated  based  on  the  PC-ADP  backscatter.  The  red  rectangle  marks  the  approximate 
sediment-flux  convergence  moment.  During  this  period,  a  ID  vertical  bed-water  column  sediment 

balance  model  is  likely  invalid. 
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